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OH  THE  PROPAGATION  OF  SMOKE  FROM  FACTORY  STACKS 


by 

L.  S.  Gaadia  and  R.  E.  Soloveichik 

I.  A  theoretical  analysis  of  smoke  propagation  froa  factory  stacks 
aakes  it  possible  to  predict  the  degree  of  air  pollution  by  smoke  in  an 
area  with  specific  characteristics  of  the  stacks  and  smoke,  which  is  ejected 
into  the  atmosphere  and  under  certain  meteorological  conditions.  By  knowing 
the  latter,  it  is  possible  on  the  basis  of  a  theoretical  analysis  to  deter¬ 
mine  such  characteristics  of  stacks  and  smoke  with  which  the  smoke  concen¬ 
tration  does  not  exceed  the  limits  of  certain  established  norms. 

In  an  analysis  of  smoke  diffusion  it  is  important  to  take  into  consider¬ 
ation,  generally  speaking,  that  smoke  particles  are  not  suspended  but 
possess  a  characteristic  (proper)  velocity.  Frequently,  as  a  result  of  over¬ 
heating  the  smoke  will  possess  an  ascending  velocity  at  a  certain  initial 
phase  of  its  propagation.  However,  such  an  initial  phase  is  of  a  short 
duration  inasmuch  as  under  the  influence  of  the  turbulent  beat  conductivity 
temperatures  of  the  smoke  particles  and  the  surrounding  air  rapidly  become 
equal.  It  is  extremely  difficult  to  calculate  accurately  the  described 
effect  because  it  would  be  necessary  to  solve  simultaneously  the  diffusion 
equation  and  the  equation  for  the  free  convection  of  smoke  particles. 

However,  even  with  a  considerable  overheating  of  smoke  this  effect  may  be 
calculated  approximately  by  substituting  the  actual  source  of  smoke  with  a 
fictitious  -  a  somewhat  elevated-source »  (see  C^jf)  *  Most  frequently 
the  effect  is  negligibly  small  due  to  the  short  duration  of  the  mentioned 
first  phase  of  smoke  propagation,  and  at  times  it  is  totally  absent.  There¬ 
fore,  vertical  velocities  caused  by  the  overheating  will  not  be  taken  into 
consideration.  The  descending  characteristic  velocities  of  smoke  particles 
which  are  caused  by  the  weight  of  particles  are  much  more  essential.  For 
particles  of  a  certain  size  the  rate  of  their  fall  w  may  be  assiaed  as 
constant.  Although  the  usual  values  of  v  are  extremely  snail,  their 
influence  nay  be  significant,  as  long  as  it  is  realized  throughout  the 
entire  tine  the  particles  are  in  notion  with  the  exception,  perhaps,  of  the 
mentioned  initial  phase.  In  addition,  the  existence  of  the  rate  of  fall 
v  leads  to  certain  qualitatively  new  effects. 

Thus,  if  the  smoke  admixture  is  polydlspersed,  i.e.  the  smoke  particles 
emerging  from  the  chimney  are  of  different  sizes,  then  during  the  process 
of  propagation  of  an  admixture  there  occurs  a  re-distribution  of  particles 
according  to  the  size  -  the  heavier  particles  fall  at  a  distance  closer  to 
the  source  than  the  lighter  ones.  Obviously,  this  fact  would  he  impossible 
to  describe,  if  the  characteristic  rate  of  fall  of  particles  is  ignored. 

The  presence  of  proper  velocities  v  leads  also  to  another  effect 
which  at  the  first  glance  is  paradoxical:  with  an  increase  of  the  strength 
of  the  source  Q,  i.e.  the  amount  of  smoke  entering  into  the  atmosphere 
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per  unit  of  tine,  smoke  concentration  over  a  given  distance  from  the  source 
increases  only  19  to  a  certain  Unit  but  vith  a  further  increment  of  Q, 
diminishes.  This  paradox  is  simple  to  explain,  if  it  is  taken  into  account 
that  vith  tbs  lac  resent  of  Q  the  average  sise  of  snobs  particles  should 
also  increase,  and  consequently  also  the  naan  value  of  v«  The  inc resent 
of  Q  alone  leads  to  the  increase  of  the  concentration  of  adnlxture  q  at 
any  distance  frosi  the  source,  and  the  laereaae  of  v  -  leada  to  the 
decrease  of  q.  At  snail  values  of  Q  the  direct  influence  of  the  in- 
c resent  of  Q  is  sore  substantial,  so  that  q  increases.  However,  with 
largo  Q  values,  the  influence  of  the  corresponding  Increase  of  the  proper 
rate  of  fall  v  is  sore  substantial  and  therefore  with  the  increase  of  Q 
the  concentration  of  q  decreases. 

The  theory  of  the  propagation  of  snoke  adnlxture  in  the  atmosphere  and 
of  analogous  processes  van  event ned  in  a  number  of  works.  Sutton  7 
and  M.  B.  Berliand  C'lJ  have  not  taken  into  account  the  characteristic 
rate  of  fall  of  snoke  particles  (see  also  /~3 J)  .  Osanedy  £\J  proposed 
to  take  into  account  approximately  the  individual  velocity  of  v  by  the 
■sans  of  substituting  the  axis  of  the  source  (l.e.  vith  a  horiaontal 
straight  line  vhich  is  directed  fraa  the  source  vith  the  vied),  vith  a 
straight  Una  inclined  towards  the  sods  downward  at  a  certain  angle.  Bvi- 
dently,  such  e  asthod  la  ne  longer  valid  already  awing  to  the  fact  that  the 
line  of  the  ■» values  of  q  differs  substantially  fraa  the  straight 
lina  in  case  vjfc). 

N.  1.  IOdla &J,  f6j,  as  veil  as  L.  B.  Bandln  and  A.  S.  Dubov  pfj 
have  exanlned  a  problem,  vhich  essentially  is  equivalent  to  the  problem  an 
two  dimensional  propagation  of  heavy  admixtures,  by  icsorimg  the  tertsontml 
turbulent  exchaiqe  or  the  sans  as  the  problem  of  propagation  of  admixture 
tram  a  linear  source .  A  similar  problem,  ves  examined  by  FOrtak  £a J  but 
under  a  more  general  initial  condition,  vltb  regard  to  the  problem  am  the 
propegatlom  of  sand  fraa  the  desert  into  the  sea,  vhareupon  contrary  to  the 
authors  mentioned  earlier,  Fortak  United  the  problem  to  the  aijq>lest  as- 
snqrtiom  an  the  constancy  of  the  vertical  turbulent  coefficient  exchange 
vith  height.  The  space  problem  on  the  propagation _of  heavy  admixtures  tee 
beam  —a  by  Booaaquet  and  otters  £9l7»  L^J*  ,rt>°  however,  did  not 
obtain  a  ccnplste  solution  to  this  problem.  Ibis  complete  solution  has 
been  found  recently  by  A.  I.  Denisov  who,  by  the  vay,  established 

the  error  in  one  of  the  Bosanqust's  formulas. 


The  present  work  is  dedicated  to  the  solution  of  the  problem  vhich  la 
concerned  vith  a  series  of  acre  general  relationships  and  vhich  descrlbad 
more  closely  the  behaviour  of  nscke  under  atmospheric  conditions  then  the 
previously  solved  problems.  In  the  analysis  of  the  obtained  eolation 
attention  is  devoted  to  specific  results  caused  by  the  effect  of  the  proper 
rates  of  foil  of  snoke  particles. 


I 


*!»*•  i  Hi  y  mtki  borlsontal  coordinates,  with  *  coinciding  with 
tha  wind  direction  Mi  y  perpendicular  to  tha  wind  respectively,  %  -  la 
tha  haltfrt,  by  -  is  ths  coefficient  of  tha  horiscntal  turbulent  diffusion, 
ha  -  is  ths  coefficient  of  turbulent  diffusion  along  the  vertical,  u  -  is 

the  wind  speed. 

I*  d  -  is  ths  concentration  of  particles  with  a  certain  characteristic 
sine  (radius)  r,  then  ths  characteristic  rate  of  fall,  of  w  ,  my  be 
assunad  as  constant  in  agra  onset  with  what  was  said  previously.  The  wind 
velocity  u  and  the  coefficient  k~  shall  be  coaside  red  as  poser  func¬ 
tions  of  height 

H  =  kK.  =  k„Zm  (2) 

with  the  seen  exponent  of  tha  poser  a. 

Such  an  amxexlmtiaa  has  been  applied  is  a  nunher  of  works  (for 
L^Jt  L^J •  £13-7)  proved  to  be  successful.  It  should  be  noted 
that  a  la  always  a  proper  fraction.  For  practical  purposes,  as  a  rule 

0.x  <.<-0.2. 

the  coefficient  kg  will  be  Manned  to  increase  linearly  with  height 

k.  =  At,  z.  (3) 

Let  us  solve  the  problm  la  the  interval  »<-*<*•.  -  -v<_y<  foo,  o  <r<  'o. 
As  the  Kneed  ary  conditions  we  shall  Msms  tha  linitatlon  of  concantratioa 
at  infinity 


f '*~oo^«e. 

(4) 

IV 1*  ~±oo-^  **’ 

(S) 

ths  give 

a  distribution  of  concentration  in  tha  "initial"  plane 

(6) 

and  1st  es  also  assent  that  tha  turbulent  flux  of  adnirtura  an  tha  under¬ 
lying  surface  Is  equal  to  nsro 

<7) 


Lot  us  asto  that  uadar  condition  (7)  the  total  flux  of  adnixture  differs 

free  aero  owing  to  ths  oadlnantation  flux  however,  it  is  United. 

It  my  be  indicated  that  for  our  problm  condition  (7)  Is  equivalent  to 
the  o sedition  of  ths  linitetlen  of  eonesntxetlon  on  ths  underlying  surface. 

(») 

®nder  such  farm  let  lens  tbs  problm  on  ths  propagation  of  heavy  par¬ 
ticles  af  anoke  is  closest  to  the  problem  of  Bosaaqast  •  Fbarsoa,  solved 
by  Beal  sow  but  differs  free  it  in  the  following  relationships . 
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1.  Instead  of  tbs  "initial”  condition  for  the  typo  of  source,  condition 
(6)  of  tbs  general  typo  Is  ussd.  It  should  bo  noted  that  the  source  of 
woks  in  tbs  fora  of  an  opening  of  a  nobs  stack  nay,  generally  speaking, 

be  considered  as  a  point  source,  l.e.  to  as suns  that 

F(y.  z)  =  Q*(.y)*(z-A).  (f) 

where  h  -  is  tbs  height  of  tbs  source,  g_  is  tbs  delta-function*. 

However,  for  a  amber  of  other  practical  problems,  as  for  example  for 
the  problem  mentlomed  above  concerning  the  propagation  of  sand,  it  is 
necessary  to  use  tbs  condition  in  the  initial  plane  of  a  more  general  fom 
than  (9);  therefore,  it  would  be  useful  to  obtain  a  solution  under  the 
condition  (6)  especially  as  the  method  for  the  solution  used  in 
under  the  condition  of  a  source  type  cannot  be  applied  directly  in  the 
case  of  condition  (6) .  Of  course,  from  the  solution  under  the  condition 
(6)  it  is  easy  to  obtain  a  "source"  solution,  for  which  it  is  sufficient 
to  determine  F  (y,  a)  by  equation  (9)  and  to  take  advantage  of  the  basic 
characteristic  of  the  delta-function. 

« 

J  <!»  (x)  S  (x  —  c)  dx  —  d*  (c),  a<c<<>-  (1°) 

a 

2.  The  wind  speed  u  and  the  horisontal  exchange  coefficient  k_  are 
not  considered  as  constant  but  increase  with  height  according  to  (2) .  It 
is  clear  that  this  gives  the  best  approximation  0f  tbs  actual  atnoepbarlc 
conditions. 

3.  Tbs  cosfflclent  k_  is  as  eased  to  be  independent  of  the  horisontal 
coordinate  x  ,  while  in  works  [~9_J  .  ['llj'  it  is  assunad  that  it  is 
proportional  to  x  .  The  last  assumption  used  in  a  number  of  worts  by 
the  English  authors,  is  based  oa  tbs  feet  that  the  intensity  of  tbs  turbu¬ 
lent  nixing  of  tbs  admixture  should  depend  on  the  characteristic  dlnaaslcms 
of  tbs  admixture  "cloud"  increasing  with  their  growth;  since  the  admixture 
"cloud"  expands  in  tbs  direction  x  then  they  assme  that  k_  increases 
with  an  lac  resent  of  x  .  Such  an  assumption  is  not  vnry  convincing  owing 
to  the  following  reasons. 

Firstly,  if  it  would  have  been  possible  to  agree  with  the  arguments 
presented  above,  then  they  would  pertain  not  only  to  the  coefficient  ky 
but  equally  to  the  coefficient  kg,  ,  so  long  as  tbs  admixture  "cloud” 
expands  with  the  increase  of  x  along  the  vertical  also.  In  addition  to 
this,  these  coefficients  should  be  considered  as  independent  of  the  x 
coordinates  but  dependant  oa  the  distribution  of  tbs  diffusing  adnixture 
q(x,  y,  s),  i.s.  s  non-linear  prefclen  should  bo  solved.  Secondly,  it  is 
not  quite  clear  in  what  manner  these  arguments  nay  be  used  in  case  of  the 
initial  distribution  of  a  general  type  (6),  lasteed  of  a  source  type. 


PM  *  t> h),  0<ir< 

-  f(U?.o  outftld*  of  tfela  lateral.  (Trau  later) 
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ttaalljr,  thirdly,  matron*  theoretical  investigation*  exist  in  which  the 
coefficient  of  the  horizontal  nixing  was  as  sued  a*  independent  of  the 
horisontel  coordinate*  (for  example  flj,  LlJ,  L^J> '  A  too* 
nent  between  the  theoretical  results  and  the  observational  data  shows  that 
the  coefficients  of  turbulent  diffusion  nay  be  considered  Independent  of 
the  horizontal  coordinates. 


m.  By  substititing  (2)  and  (3)  into  (l)  we  shall  obtain 


u  !?£  — di  _  u  .»  i  1,  dq 

u'Tx~wz  Tz  =  *>* 


Let  us  assuns  at  first  that  the  distribitlon  of  <1  in  the  initial  plane 
synaetrical  with  respect  to  the  axis  s,  i.e.,  that  function  F  (y,  t)  - 
is  an  even  function  according  to  y  .  Then  by  Multiplying  equation  (11)  by 

y  (y;  X)  =  cos ).  y  (12) 


and  as  swing  that 


J  V(Jf.  y.  -)  Y(y\  — «|-( .v.  r;  >.) 

—  O' 


(13) 


*"*  //•(>•.  r)K(y;  /)rfv-=/(z;  >.).  (M) 

—  Of, 

We  obtain  according  to  (12)  and  (13) 

JO 

q(x,  y,  2)  =.  f'Y(x,  ).)  Y( y;  i)rfX.  (15) 

II 

By  integrating  equation  (ll),  which  is  Multiplied  by  Y  ,  fix*  -  ^  to  <x, 
with  respect  to  the  variable  y  ,  we  shall  obtain,  by  using  condition  (5), 

..  _ «.  -m  d  /_  <W\  .  i*  0‘Q  . 


ll  —  °'1  b  >  “  m  **  u  x,  A 

“•  dj - wz  -dt = a (* Tz) ~ *» ./  ^  y dy 


or  after  integrating  twice  by  parts  in  the  last  teza 


dx~WZ  dz\Z'Oz) 


Further  let  us  perfom  an  operational  transfomatlon  according  to  x  . 
Since  according  to  (6) ,  (13)  and  (l4) 

*F  I  =  f . 

then  the  transfonation  of  equation  (16)  gives 

">P 0 *  -  f)  -  wz-  f^  kxz  mdJz%)  ~  W*. 

or 

tf-M  ,  /.  ,  a>\  1  */'•  *11*-  '  u, p  m  1,7.  _«i 

d* + 1 1  +*r/7  rf*  2  -  *, z  J ' 


Where  f  —  is  the  transfom  of  function  *1' 


* 


The  general  solution  of  equation  (l?)  nay  be  written  out  In  the  fora: 
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where  A  and  B  -  are  constants  of  Integration.  For  determining  then  we 
shall  use  the  boundary  conditions  (4)  and  (8),  which  lead  to  the  conditions 
for  qT 

As  a  result  of  the  first  of  these  conditions  A  »  0,  and  owing  to  the 
second  B  *  0,  therefore  the  solution  for  V  obtains  the  font 


Owing  to  (l8) ,  (k\z^)z  o  =°  •  therefore,  the  solution  also 

satisfies  condition  (7)  • 


In  order  to  obtain  the  original  '*r  (*•  71  fron  the  transfozn  V  (z;  p,  *) 
let  us  use  the  well  known  formula  /"see  for  instance  /14_/>  formula 

(9.  140)  J.  _  _>±1 


where  t  in  this  case  is  the  coordinate  x. 
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X&. 1  **  If-  v  le'-* 


by  canblalag  the  integrals  w  shall  obtaia 

*-it+%I5r  '  //«: ‘K'ltf  "  """  x 

II 

Beturalag  to  q  ,  n  Obtain,  according  to  (12) ,  (it) ,  (15)  aad  (19) 

OO  OO  OO  hn\*x  w 


0  -oou 

u,(;i 

x«“  (,+ “’**■'  /  . 

or  integrating  with  raapaet  to  X,  finally 


I  -fw 


2Bi(*C) 


(1  +  m)*ktx 


di  dt\  dr. 


4(l+*)/«M1  {*)  f  J :) 


(]  — oc 


_  Mtiy-tf 

e  Ujr  + 


_  «.(y+m> 
+e  'K* 


9  m,  (t1  +«j-  ,*  +  *) 

!+•  | 

C-fiV*'*  o+^M  7 

2u,(*C)  2  I 

(!+•)*, 

(1  +m)2  k,x 

The  probloa  would  ba  solved  analogously  vbaa  7  (y,  a)  is  tha  odd 
function  of  tha  eoordlaata  y  .  Tbea  sinky,  should  be  takea  as  a 
factor  for  X  aad  according  to  this  tha  final  formal*  will  differ  fras 
formula  (20)  only  by  tha  fact  that  Instead  of  the  stsa  of  tha  exponent  tha 
difference  will  appear  la  square  brackets,  flaally,  la  gsaeral  case  it 
Is  staple  to  obtaia  the  solution,  by  taking  tha  function  F  (y,  s)  ia  tha 
fox*  of  tha  son  of  the  area  aad  tha  odd  function 

F{y.  *)=  \  I F{y,  z)  +  F(—y>  *)I  4-  7 \F(y<  *)-F(-y,  *) I- 


wa  obtaia 

9  = 


= — v=j  f  7  jW*-v+ 

4(l  +  <rlt/**0*l  *'  J  J 

<1  —jo  L 


v  U  +  Ml l«M  /  I  *»i  <f -I _ 

'  •  It  I  + 

n  *1  *1  1 


dtjd:. 


(21) 
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Formula  (2l)  gives  a  complete  solution  of  the  stated  problem.  In 
particular,  for  the  point  source  (9)  ve  obtain  from  (2l)  or  (20)  on  the 
basis  of  (10) 


q 


2  ( l  +  m)  y/nt,  *, 


«.>* 

4  kji 


u,  ~  m) 

H  **•*,» 


/  . 


1  *■  m 

■la,  (A*)  1 

(T+«VM  ■ 

i 


(22) 


From  (22)  are  obtained  simple  expressions  for  the  distribution  of  the 
axial  concentration  ^  la  case  the  source  is  elevated  (y  =  o,  z  =  h) 


V*  = 


2  <H-/»)V/i»*u*1  8 


2u,Hl  +  m 
ll-  m?k,z 


(i  —  m\ It, 


\  2uy+>"  ) 

ill  t  m)‘k,x  j 


(23) 


and  for  the  distribution  of  the  axial  concentration  q©  in  case  the 
source  Is  at  ground  level  (y»2  =  o> 


„,n 1 


<?o=  - 


(*  QH  *' _ (■ _ u, _ -,i 

'"»r(i+(n%*1l^*-  (1 


(24) 


IV  In  analyzing  the  solution  it  is  convenient  to  change  to  dimension* 
less  variables  by  means  of  relations 


1  =  ti _ _ =  „  _ 1  /  ^  h'~m  q 

V  ky<c  J  ^  (Itw)  k,  P‘  1/  *■  " 


_X_  i/**t 

(l»m)  |  *1 

1  - 


4sVh’”"-~ 

Then  instead  of  formula  (22)  we  obtain 


«(t,  C:  //)-//•("/* 


-  //»- ; 


(23) 


(20) 


and  correspondingly  instead  of  formula  (23)  and  (24) 

s(0,  //;  W)  -.«•*““•/,  (2tf») 

•(0,  0;  AO-rTr|p7y«'+V«' 


(27) 

(2*> 


Formula  (27)  shams  that  the  axial  concent  ration  for  am  elevated  source 
decreases  aonotonlcally  with  an  increase  in  the  distance  from  the  axis 
of  the  source,  l.e.,  with  the  deerease  of  B  .  Along  all  the  other  hori- 
» octal  straight  lines  in  the  plane  y»0  the  adr’xture  concentration  at 
first  increases  with  distance  from  the  point  x-0  ,  and  only  after  reeoh> 
lag  a  certain  earl  mm  begins  to  decrease.  At  the  axis  of  the  ground 
sourco  is  particular,  the  maximal  vales  «,  according  to  (20) ,  is  rose  bad 
at  the  point 

(29> 


I 


and  is  equal  to 


o.=o(0,  0;  HJ 


l 

l  0+/» 
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(30) 


It  is  interesting  to  note  that  the  distance 

_ 1 _ _«i 

r  — .  i  ' 

(i  +m)- 


-1 


(i  +*j*, 
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If 


(31) 


at  vhich,  according  to  (29)  and  (25)  the  waxlnun  concentration  is  reached, 
is  independent  of  the  horizontal  diffusion  coefficient  Isq  ,  and  the 
naxinun  concentration,  being  equal,  according  to  (30)  and  (25) 


(I  4  my 
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0/*, 


(1  + 


_+l 

ll».  1 


(32) 


is  independent  of  the  wind  speed  . 

Both  of  these  conclusions  are  clear  from  simple  physical  considerations. 
Indeed,  the  naxlnnn  concentration  should  be  dependent  only  on  tine  JL 

ui  ’ 

in  which  particles  are  transported  by  the  vlnd  fron  the  source  to  the 
point  of  the  abscissa  x  ,  and  this  tine,  apparently,  is  not  dependent 
an  ui  .  Regarding  the  coefficient  k*,,  its  value  should  effect  the  rate 
of  naxinun  concentration,  but  at  the  point  location  where  such  a  concen¬ 
tration  is  reached. 

Let  us  note  that  the  presence  of  a  point  (source)  with  a  naxinun  con¬ 
centration  at  the  underlying  surface  is  not  characteristic  for  the  dif¬ 
fusion  of  heavy  particles  only.  In  the  case  of  weighted  particles  the 
general  character  of  the  distribution  of  q  on  a  straight  line  y*x»0 
is  the  sane.  In  this  respect  the  influence  of  the  proper  weight  of 
particles  of  the  admixture  is  of  a  quantitative  nature  only.  Kanely,  as 
nay  be  seen  from  (31)  and  (32),  with  the  increase  of  the  deposition  rate 
of  particles,  the  point  at  vhich  the  naxinun  concentration  is  reached, 
approaches  the  source  but  the  value  of  the  naxinun  concentration  increases. 

In  general,  the  special  distribution  of  a  concentration  does  not  undergo 
quantitative  changes  under  the  Influence  of  the  weight  of  admixture 
particles.  Therefore,  we  shall  not  describe  the  distribution  of  q  any 
further  but  we  shall  turn  to  the  analysis  of  specific  effects  caused  by  the 
influence  of  the  proper  rates  of  fall  of  particles. 

Let  us  assume  that  the  admixture  is  polydlspersed  and  the  distribution 
curve  for  the  number  of  particle 8  emitted  by  the  source,  according  to 
the  sixes  of  r  is  described  by  the  function  f (r) .  Tor  simplification 


we  shall  assume  the  particles  to  be  spherical  and  r  vill  be  considered 
as  the  radius  of  the  particles.  These  considerations  nay  easily  be 
generalized  in  the  case  of  particles  which  are  arbitrary  in  shape,  pro¬ 
viding  there  is  no  relationship  between  the  size  and  the  shape  of  particles. 

If  I0  is  the  total  number  of  admixture  particles  emitted  by  the 
source  per  unit  of  tine,  then  according  to  what  has  been  said, 

Ntf(r)  dr 

is  the  nuaber  of  particles  (among  the  total  nuaber  HQ)  with  a  radius 
between  r  and  r  +  dr  .  Then  the  strength  of  the  source  Q  is  equal, 
obviously, to 

Q  =  /  4«W(r)rfr. 

where  p  —  is  the  density  of  admixture  substance  and  the  integration  is  made 
throughout  the  entire  spectrum  of  the  size  of  particles. 

Let  us  assume  that  the  size  of  every  particle  does  not  change  during 
the  diffusion  process  of  the  admixture  and  that  no  interaction  occurs 
between  the  particles.  This  asststption  nay  lead  to  certain  errors  only 
at  the  1  Mediate  proximity  of  the  source,  but  at  a  distance  from  it,  where 
the  concentration,  q  ,  is  small,  it  is  very  close  to  reality.  Following 
this  assumption  it  is  possible  to  examine  the  behaviour  of  particles  of 
every  fixed  size  independently. 

Let  us  substitute  in  (28) 

3(0,  0,  H)  =  ^-S, 

in  order  to  single  out  the  dependency  of  the  concentration  on  3. 

Then  formula  (28)  for  the  concentration  of  particles  of  a  given  size 
at  the  ground  level  axis  of  the  source  vill  be  written  in  the  form 

where,  in  addition  to 

Q,- ^.«rVV,/(r)  .  (34) 


also  only  the  parameter  p  ■  pr  depends  on  the  radius  r  of  the  particles, 
so  long  as  it  is  related  according  to  (25)  to  the  deposition  rate  of 
particles,  lamely,  owing  to  the  mealiness  of  particle  slses  it  is  natural 
to  assume  their  deposition  rata  proportional  to  the  square  of  the  radius 


P= 


<*’0 

(I  +m  )*,rj 


then 


(35) 


r 


t' 


By  having  a  eflKnte  distribution  ctm  f  (r)  of  particles  which  are 

emitted  by  the  source,  according  to  sixes  end  knowing  the  parameter  ^ '  , 

0 

It  ie  possible  on  the  basis  of  (33),  (3b)  end  (35)  to  find  the  distribu¬ 
tion  of  the  concentration  of  particles  of  given  sises  at  the  underlying 
surface.  After  this,  integrating  with  respect  to  r  ,  it  is  also  possible 
to  find  the  distribution  of  total  concentration.  In  addition  to  this,  it 
is  essential  to  keep  in  niad  that  the  distribution  of  total  concentration 
will  vary  from  the  distribution  of  the  corresponding  nonodispersed  ad¬ 
mixture,  which  is  obtained  by  substituting  the  sises  of  all  particles 
with  the  average  sites.  In  particular,  at  the  earth's  surface  the  total 
concentration  of  polydispersed  admixture  mil  be  greater  than  the  cor¬ 
responding  concentration  of  the  nonodispersed  adalxture. 

It  is  even  note  essential  to  take  into  account  the  re-distribution  of 
particles  according  to  the  else  in  the  propagation  process  of  adadxtnre. 
By  considering  the  concentration  at  the  underlying  surface,  it  is  natural 
to  introduce  the  coefficient 

H^r 

*<r)~r(l  '  (36> 

which  describes,  according  to  (33)  the  distribution  of  s  according  to 

40 

else  if  -jp  —  0,  i.e.,  in  the  event,  the  strength  of  the  source  for 

various  sixe  particles  is  the  sane.  In  other  words,  the  coefficient  x(r) 
describes  the  influence  of  the  deposition  of  particles  on  the  distribu¬ 
tion  of  their  mass  according  to  their  sises  over  the  underlying  surface. 

The  behaviour  of  function  x(r)  differs  depending  on  H  .  Handy,  for 
smell,  ■  values  (greater  distances  of  x  )  /  (r)  has  a  monotonlc  decrease 
with  the  increase  of  r  ,  but  for  large  H  values  (small  x  values)  x  (r) 
increases  to  a  certain  swzlnun  and  only  after  that  it  decreases.  It  is 
simple  to  describe  this  effect  quantitatively,  if  we  keep  in  aind  that 

/>,<<!  (due  to  the  smallness  of  the  parameter  |TnTmT*T]  )  •  The  con¬ 
dition  of  the  naxlnua  of  z(r)  gives 

2  In  H  —  <|»(i  p)  «*  0, 

where  <j>  —  is  the  logarithmic  derivative  of  a  gamut-function. 


1  If,  for  example,  the  deposition  of 
Stokes  law,  then  w 


particles  takes  place  according  to 

9h 
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By  the  expensive  of  the  function  lato  a  oerl* 
two  tons 

•+  6P  ’ 

wo  obtain  the  value  of  p  ,  corresponding  to  the  ■ 


Since  always 
only  if 


0.-.T  ('■«+!). 

P>0,  then,  according  to  (3?)  tbs 


by  confining  to 

of  x«- 

(37) 

■  of  X(r)  occurs 


i.e.,  in 


with  (25),  if 

*<•*» 
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(t  +m  »- 


(38) 


A  physical  explanation  of  this  fact  consists  in  the  following.  Gener¬ 
ally  speaking,  the  larger  are  the  particles  than  to  a  greater  degree  is 
their  concentration  decreased  at  a  given  distance  frea  the  source,  because 
a  large  portion  of  such  particles,  owing  to  deposition,  fall  oat  at  closer 
distances.  But  at  short  distances  (38),  the  concentration  of  very  fine 
particles  is  anall  owing  to  the  fact  that  as  a  result  of  the  slew  depo¬ 
sition  rate  they  do  not  quite  reach  the  underlying  surface.  At  noderate 
distances  with  all  conditions  being  equal,  there  ere  eostly  particles  of 
intenusdiate  aixas  which  satisfy  foraula  (37) . 


It  is  interesting  to  note  that  the  upper  Unit  xm  of  the  distances, 
where  the  curve  /  (r)  reaches  the  naxlaon,  is  very  slaply  associated  with 
the  abscissa  x  ,  for  which  tha  concentration  of  particles  of  a  osrtaln 
given  also  la  naxleal.  According  to  (31)  and  (38), 


e 


w 

(1  +  ™)*i 


^4,38. 


Finally,  let  us  sxsnlna  tha  p  rob  leu  concerning  tha  relationship  be¬ 
tween  the  concentration  of  particles  and  tha  strength  of  tha  soeroe .  Xt  la 
Obvious  that  if  the  process  of  tho  propagation  of  adulxtars  is  independent 
of  tho  strength  of  the  source,  than  at  each  point  tho  concentration  la 
proportional  to  the  strength  of  the  eonree  end  increases  without  Halt 
with  the  increase  of  the  strength  of  the  source.  Actually,  however,  this 
la  aors  complicated.  Sanely,  with  the  chsngo  of  the  strength  of  tho 
source,  the  distribution  curve  of  the  else  particles  being  enitted  by  it 
also  changes,  so  that  with  an  lacraaea  in  the  strength  of  tho  soeroe,  the 
average  particle  sloes  Increase  also.  Therefore,  the  yaxtlelee  settle 
■ore  intensely  and  it  can  ba  asswnsd  that  at  a  certain  far  distance  trm 
the  source,  as  a  result,  the  concentration  with  the  increase  ef  the 
strength  of  the  source  dees  net  grew  but  dial al shoe. 


of  such  bo  effect  end 


Let  us  enod.Be  the  slnplsst  physical  schwas 
aaaely  lot  us  as  stew  that 

a)  the  adaixture  Is  aonodlspersed  sad 

b)  the  amber  of  particles  which  is  ealtted  by  the  source  per  unit 
of  tine  is  independent  of  the  strength  of  the  source. 


Then  it  is  obvious  that 


where  Qq  -  is  the  value  of  Q  at  r«r0. 

Owing  to  (35)  sad  (39) 

.3 

Q=Qo('ji  IK-Pr-,) 

By  substituting  (to)  into  (33)  we  obtain 


S  —  _H3+  ip  1  -  w 

.?  nr T7TP  e 

Pi 
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The  condition  of  the 
cording  to  (tl) , 


ustI—  of  s  with  respect  to 

3  -i-  ♦/>  In  //  —  2p 7  ( 1  +J>)  =  0. 


P 


gives,  ac- 


fren  which  by  confining  to  the  first  tern  of  the  expansion  of  y  (l  p) 

into  a  series 


P  =  — 
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In  H  + 


(42, 


Fran  (t2)  it  is  evident  that  s  as  a  function  of  p 
a  function  of  Q  )  has  a  naxinm  only  IF 


(or  the  sane  as 


l.e.,  if 


In  H<-{, 


-  P  u]  “T  m 

x>x«—  (i  -;-mH  *1  " 

The  lower  lisdt  •*„  of  these  x  values  coincides  with  the  upper  Unit 
of  fomula  (38).  Thus,  for  all  sufficiently  large  distances  of 
the  concentration  with  an  increase  in  the  strength  of  the  source  does  not 
increase  infinitely  but  reaches  a  certain  narlnm  with  the  Q  values, 
which  is  obtained  by  substituting  (t2)  into  (to),  and  with  a  further 
increase  Q  decreases.  The  values  of  p  and  correspondingly  of  Q  at 
which  the  indicated  anrlum  is  reached,  essentially  depend  on  the  distance 
x  fron  the  source,  and  sanely 
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